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Abstract 

We discuss the approach to equihbrium of systems governed by the 
Fokker-Planck equation. In particular, we focus on problems involv- 
ing barrier penetration and the associated Kramers' time. We also 
describe the connection between stochastic processes and quantum 
mechanics. 
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1 Introduction 

In this article, we will discuss the approach to equilibrium of a system with 
fluctuations. The simplest example of this is an overdamped particle with co- 
ordinate x{t) at time t. This particle experiences an external potential U{x). 
The time evolution of the position variable obeys the Langevin equation, 
which is a stochastic differential equation: 

where F is the inverse friction constant and f{t) is the random or fluctuating 
force. The distribution of fit) is Gaussion, and its correlation function obeys 

<fiti)f{h)>=2e5{h-h) (2) 
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where e measures the amphtude of the noise term. 

The equihbrium distribution for x will be proportional to exp[— t/(a;)/e] 
if the fluctuation-dissipation theorem holds, flxing F = e. It is well known 
that if P{x, t) is the time dependent probability distribution for the random 
process described as x{t), then P{x,t) satisfies the Fokker-Planck equation 
(Risken, 1984): 



dP{x,t) 
di 




(3) 



It is clear that if P = Pq oc exp(— t//e), then ^ = and this corresponds 
to the equilibrium distribution. The approach to equilibrium from a non 
equilibrium state will be the subject of the ensuing sections. 



2 Kramers' Time 

In this section, we recall the approach to equilibrium in a system, whose time- 
dependent probability distribution is governed by a Fokker-Planck equation. 
In particular, we discuss the situation where the potential driving the dy- 
namics is bistable. For the most part, we will talk about a one-dimensional 
potential. 

A convenient way of handling the problem is to make the transformation 

P(x, t) = exp(-f//2e)(/)(a;, t) (4) 
For the new variable 4>{x,t), one obtains the Schrodinger-like equation 

where 

1 U'ix')'^ 1 

V{x) = - tu"ix) (6) 

If (pnix) are a complete set of eigenfunctions of the Hamiltonian H — 
— + then we can expand as follows: 

t) = ^an exp{-Xnt)(j)n{x) (7) 

n 

where the A„ are the eigenvalues of H given by 

H4>n = K4>n (8) 

Further, the constants determined dy 0), which is obtained from 

P{x, t). It is clear from the form of H, which can be written as 

H = eA^A (9) 
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with 



Figure 1: The general bistable potential in one dimension. 



A- — — 

dx 2e 



(10) 



that the spectrum of H is non negative and the lowest eigenvalue is zero 
with the corresponding eigenfunction given by 0o = Nexp{—U/2e), where 
N is a normalization constant. It is clear from Eq.((Zj) that as t — > oo, only 
the term n = will servive (Aq = 0) and the limiting value of (j){x, t) is 
exp(— f//2e), corresponding to equilibrium distribution Peq{x) oc exp(— f//e), 
as is obvious from Eq.Q. Thus, in the Fokker-Planck equation, the approach 
to equilibrium is guaranteed. 

The approach to equilibrium occures on a fast time-scale if the potential 
has a single minimum, e.g. U{x) = x"^. In this case, V{x) = — 1, A„ = 2n 
(n = 0, 1, 2, ) and Eq.Q reads as 



0(x, t) = ao(j)o{x) + ai0i(x) exp(— 2t) + 



(11) 



The exponential term dies out on time-scales of 0{1) and the system, 
if started out from a non-equilibrium state, will evolve to the equilibrium 
state on a fast time-scale. Specifically, if we take the initial probability 
distribution to be Po{x) = exp(— (x — a)^/e), then 0(a;, t = 0) = 

l/(7re)"'^/^ exp(— x^/2e + ax/e — a^/2e), with the eigenfunction 0„(a;) given 
by 0„(x) = {-T/ 



^172^172^ )^^^exp(-xV2e)iy„(x/e^/2)_ ^-^^ expansion 

coefficients a„ of Eq.((Zj) to be given by a„ = io,/^^^'^)^2^. 

Eq. (IZj) is easily performed keeping in mind the generating function of Hermite 

polynomials to yield 



P{x,t) 



1 



(vre)!/' 



exp 



[X 



ae 



-2t\2 



(12) 



We now turn our attention to the bistable potential shown in Fig.(l). At 
t = 0, we choose a probability distribution centered sharply around x = 0, 
which is clearly a non-equilibrium function. Because of the inverted oscillator 
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potential near the center, the distribution will broaden out in the initial 
stages. Thereafter, in an intermidiate time zone called the Suzuki regime, 
the broadened central peak begins to split into two side peaks corresponding 
to the minima of U{x) at x = a and x = b. Finally the system enters the 
Kramers' (1940) regime where it is close to equilibrium but makes occasional 
large fluctuation due to the noise. This Kramers regime will be governed 
by the lowest eigenvalue Ai of the system. The inverse of Ai sets the time 
scale for attaining equilibrium. The eigenvalue Ai is exponentially close to 
the ground state value Aq = and differs from it due to the tunneling in the 
three-well problem. A WKB calculation of the eigenvalue was carried out by 
Caroh et al. (1979). 

Here we will show a variational calculation (Bhattacharjee and Banerjee, 
1989) motivated by the pioneering work of Bernstein and Brown (1984). The 
observation of Bernstein and Brown was that the Hamiltonian H = eA^ A and 
H = tAJ^ are surpersymmetric partners having the same spectrum except 
for the ground state. While the ground state energy of H is zero, the same 
can not be true for H because the wave function satisfying = cannot 
be normalized. The ground state of eAA^ is then the first excited state of 
eA^A. This is easily seen if we construct the matrix 



The operator Q with the representation ( ^ !! ) commutes with H. 



s 



A 

and hence Hg has a degenerate spectrum. If ip is an eigenfuntion of Hg with 
eigenvalue A, then so is Qip. By inspection, one set of eigenvalues of Hg can 



be written as ipn = ( g" ) since (see Eqs.® and 



H. ( t ) = ^ ( ■'*0*' ) = ^» ( t 1 (") 



Consequently, we must have 

Since ( I = f . . , I , it follows that 

V A(f)„, J \ eAA^AcPn J 

eAA^A(j)n = XnA(j)n (16) 

Thus, A„ are the eigenvalues of eAA'^ with eigenfunctions This will 
be true for all 0„ except 0o- The eigenvalue zero is nondegenerate. For all 

other 0„ we have ^ ^ and ^ '^"^ ^ as the degenerate eigenfunctions of 

Hg and write 0„ are the eigenfunctions of eA'^A; the A0„ are the eigenfunc- 
tions of eAA^ . To find the first excited state of eA^A (the one corresponding 
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eigenfunction we need to find the ground state of eAA^ . Now tliat tlie 
problem is one of finding a ground state energy, tliis can be easily done by 
variational calculations. 

We begin with the observation that the nonnormalizable wavefunction 

exp ^f/(a;)/2ej satisfies eAA^exp {u{x)/2e^ = 0. This motivates the choice 
of the wavefunction in the form exp (ip{x)/2e] with 



U{x), for 6(1 - /3ei/2) < x < a(l - 

U{a - aae^^"^) - U'{a - aae^/^)(x - a + aae^/^), 

forx>a(l-aei/2) (17) 

U{a - a(3e^'^) + U\a - a(3e^/^){x - b + l3e'/^), 



In the above, we have assumed e <^ 1 and anticipated that the distance 
scale to beO(e^/^), so that the variational parameters a and /? are numbers 
of 0{1). Note the matching occurs near the turning points in the classical 
regions for the corresponding Schrodinger equation. We need to evaluate 

J —oo 

The calculation involves the following steps : 

(i) Evaluating the normalization integral. We note that ip{x) ~ f/(0) — 
^U"{0)x^ for X < 1 and the maximum value of ip{x) dominates the integral 
for e ^ 1. Thus, the normalization integral is e^°^^{2TTe/U"{0)). 

(ii) The kinetic and potential energy term completely cancel iplx) = U{x). 
Thus the integration in the numerator of Eq. p8|) involves integrating from 
a(l — ae^/^) to oo and from 6(1 — /5e^/^) to — oo. 

(iii) In the range of integration discussed in (ii), it is sufficient to approx- 
imate the potential U{x) by a quadratic expression for e <^ 1. Terms like 
U{a — OLa^l'^) and U'{a — aae^/^) need to be expanded to 0(e), keeping in 
mind U'{a) = U'{h) = 0. We thus obtain \i{a,p). 

Minimizing Ai with respect to a and f3 leads to the conditions: 

a^a^U"{a) = b^p^U"{b) = 3 (19) 

Thus, we find 

Ai = c[| f/"(0) I /f/"(a)'/'e"^"/^ + | f/"(0) | /U"{bf^^e-^'/'] (20) 
with Aa^b = U{0) - U{a, b) and 

1 6^/2 

For the WKB result of Caroh et al., c = l/27r ~ 0.159. 
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Figure 2: The general bistable potential in two dimensions. 



To end this section, we consider generahzation of the above treatment 
(Bhattacharjee and Banerjee, 1989) to higher dimensions. In two dimensions, 
the bistable potential U{x,y) has the form shown in fig. (2). The important 
assumption about U (x, y) is that there exists a most probable escape path 
(the instanton trajectory) connecting the two wells. In this case, we can 
always choose the axes such that the desired path is the x-axis. About this 
path, U{x, y) can be expanded as 

U{x,y) = U{x) + ]^W{x)y'' (22) 

In the above U{x) is the one dimensional potential that we have already 
discussed. The Hamiltonian, corresponding to Eq.® can now be written as 

H = -e\^ + ^\+V{x,y) (23) 



where 

V{x,y) = ^ 



dx"^ dy'^ 



dx I \ dy 



Q^2 Qy2 



The ground state with zero eigenvalue is exp \^ — U{x, y)/2ej . As before, 
our interest is in the first excited state. 

We begin by noting that the first excited state will involve changes, mainly 
in the x-direction, retaining the y-dependence of the wavefunction as exp ^ — 

W{x)y'^ /2e^ to a good approximation. The "stiffness" in the y-direction 

suggests a variant of the Born-Oppenheimer approximation. We solve the 
y-dependent part of the above Hamiltonian, treating the x-dependence as a 
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parameter. The resulting eigenvalue will be a function of x and can be treated 
as an effective potential for the one-dimentional problem in x. Introducing 
Eq.(j22I) into Eq. lj^ . the y-dependent Hamiltonian can be written as 

Hy = -e^ + ^ _ p^//^)^2 (25) 

The lowest eigenvalue is 1{W^ + U'W - eW'fl'^ and we use this as an 
effective potential for the one dimensional problem in x. The function W{x) 
is going to be slowly varying and we can expect W'/W and W"/W to be 
small. This allows us a binomial expansion of the lowest eigenvalue and 
allows us to write the Hamiltonian as 



H 



dx'^ 



l\2 



4e 



U" U'W 



AW 



W 



d U' 

dx^Yt^m 



eW" 

w 

w 



d U' 



(26) 



where in the second step, we have deliberately not matched the higher- 
order term {W'/Wy. We require the two lowest eigenvalues of H in Eq. (j^ . 
The lowest eigenvalue is of course known to be exactly zero and so it is 
the approximate determination of Ai which is our concern and that in ac- 
cordance with our previous analysis, this is simply the ground state of the 
supersymmetric partner 



H 



d U' W 
dx^Ye^m 



d U' W 

~ dx^ Y^ m 



(27) 



We now take over the variational calculation for the one- dimensional case 
dicussed above. The only difference to be noted is that the non-normalizable 
function which gives zero eigenvalue is now 



ipi^x) = exp(f/(x)/2e) exp 



1 rw'ix] 



4 



W{x] 



dx 



(28) 



For e ^ 1, the extra term yields the factor exp(|/n-j;|7^) for the range of 

integration a(l — ae^^^) to oo and a similar term for the range 6(1 — Pe^^"^) 
to —oo, when one calculates the expectation value. This leads to the answer 



Ai = c 



U" 
U" 



1/2 



l/2-Aa/e 



+ 



EL 



1/2 



{W,/Woy/'e 



l/2-Ai/e 



where the constant c is once again 



e3/2 1 
6 (67r)i/2- 



(29) 
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To end this section we would like to point out a formal analogy (Schnei- 
der et al., 1985) between a quantum problem with potential V{x) and the 
Langevin process. In Eq.©, if we write t = ir, we get the Schrodinger 
equation (e = h/2m, V = hV) 



ih- 



+ V{x)<P 



(30) 



dr 2m dx'^ 

and the stationary states in the r-variable are characterized by the eigenval- 
ues An- In the path integral formulation of quantum mechanics, this equation 
corresponds to the measure 



(i/i(x, t) 



i/h J dr 



(m/2)(|^)2-y(x) 



~ e 



' J 2me 



iV[x{T)] 

V[xit)] 



(31) 



Turning to the Fokker-Planck equation, note that if at t = 0, a; = 0, 
then P(a;, to) = S{x — Xq) and this fixes the coffecients a„ in Eq.((7j) as = 

exp{Xnto)4>n{xo) / (po{xo) leading to 



P{x,t)=Mx)Yl 



4>o{xo) 



-A„(t— to) 



with X — Xq at t = to- The two-time correlation follows from 



(32) 



< x{t)x{to) > 



dx J dXoP{x,t \ Xo,to)xXoPeg.{xo) 

oo 

dxdxoXXo(po(x) y 



^^-^"(*-o)0„(x)0^(xo) 



n=0 



/ dxQ(\)Q{xQ)XQ(\)n[xQ) I dx(f)o{x)x(j)n{x) 

n=0 

^\<0\x\n>\^ e-An(i-io) 



For t ^ to, only the lowest eigenvalue will contribute and hence 



< x(t)x(0) >- 



|< I X I >| 



< X 1 > r e"-"! 



Alt 



(33) 



(34) 



This particular correlation function, this explores the lowest eigenvalue 
of the quantum problem. Correlation of composite variables will explore 
other eigenvalues A^. We thus arrive at the following result: If the spectrum 
of a quantum mechanical problem with potential V{x) is to be computed. 
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then it should be possible to do that by studying a Langevin process with 
potential U{x), where U{x) and V{x) are related as in Eq.®. Integration 
of the stochastic differential equation then allows one to compute correlation 
function < x(t)x{0) > as 

<x{t)x{0) >= lim ^ / x{t')x(t' + t)dt' (35) 

T->oo T J 

The important issue is to find U {x) if V{x) is given. If the ground state 
energy oi H = T + V , where T is the kinetic energy is Aq (this will not be 
zero in general), then U has to be found from 

^U\xf - {l/2)U" = V{x) - Ao (36) 
If 00 is the ground state of i7, then 

?7 = -2eln0o (37) 



3 Suzuki Regime 

In the bistable situations described in the previous section, the intermediate 
time zone where the probability distribution acquires bimodal character is 
also particularly important. We consider the evolution of an initial proba- 

2 4 

bility distribution P{x,t = 0) in the bistable potential V{x) = — ^ + ^. At 
t = 0, we assume that the probability is calculated about a; = in a manner 
which is almost like a delta function (a very narrow Gaussian in reality). 
As time evolves, the distribution is going to spread out and then develop 
two peaks. The spreading is the initial time regime, while the intermedi- 
ate time regime, where the probability distribution acquires the two peaks 
is called the Suzuki regime because of the interesting scaling observed by 
Suzuki (1977) in this time-zone. If the probality distribution explores only 
the small-x part of the available space, then in the initial time regime, we 
can approximate V{x) as V{x) = — Carrying out the transformation 
P{x,t) = (p{x,t) exp(— V"/2e),we find 




The eigenvalues of the operator H are clearly — (n + 1) with the eigen- 
function 0„ = Nn exp(— x^/e)if„(a;/2e), where Hniy) is the n*^ Hermite poly- 
nomial and Nn is the normalization constant. We can immediately write 

P(x, t) = Y, C„e-("+i)*Ar„e-^'/^if„(x/2e) (39) 

The constants C„ are fixed by the initial form of P{x,t). If we take the 
idealized form of P{x,t = 0) = S{x), then it follows that 
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1 



[27re(e2* - 1)]V2 



' 2e(e^*-l) 



(40) 



It is obvious that the probabihty distribution spreads out in time since 
the width of the distribution is given by 2e(e2* - 1) ~ 2ee^K This is the 
initial time-regime, where the distribution sees the potential as an unstable 
inverted oscillator. The validity of the result for times such that ee^* <^ 1 
(ie. the spread is not too big) or t <^ to ~ /n(l/e). 

We now consider i ^ to but much less than the Kramer's time. In 
this limit the system 'slides down' the potential and does not encounter the 
diffusion term Fokker-Planck equation. We can now write 



dP 
'dt 



Writing e = t', we have 
,dP 



d_ 

dx . 



Pi-x + x-^ 



dP 

-x + x^)— + (Sx^ - l)P 
ox 



dP 



t'— = i-x + x')— + {3x'-l)P 



dx 



(41) 



(42) 



From the method of characterisation, we can write down the solution as 



p{x,t) = ^^nh^M^vmi-x^)) 



1 ^/e 



x(l — x^) VI — X 



(43) 



where G{z) is an arbitrary function. As always G has to be obtained from 
the initial condition and in this case that is provided by Eq. (j40|) . The initial 
condition also entails x <^ 1 and comparing with Eq. ()40|) . we can obtain 



= ^(1^ Mi^tiyp -P ( - -V2(l - - 1)) (44) 

What emerges is the new time scale r = ee^*, which determines when the 
probality distribution will acquire the two peaked structure. [The sequence 
of P(r, t) for the above function with P{r, 0) ~ S{r) is shown in fig. (3).] We 
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will demonstrate that the above answer is exact if consider a A^-dimensional 
vector Xi{t), i = 1,2, ,N, with the stochastic time-evolution 



dt 



dV 

-1^— + fa 



(45) 



with V = 
correlation: 



Y + ^ = J2iLi ^1 S'lid fa is the random term with Gaussian 



< fa{t2)fp{ti) >=2eS{h-t2] 



(46) 



The fluctuation-dissipation relation holds with 7 = e. The Fokker-Planck 
equation for P{xa,t) is 



— -d (P^] c 
dt °V dxa^ dxadxc 



(47) 



For the spherically symmetric V, we expect P = P{r,t) and thus in 
spherical polar coordinate system, Eq. ()47|l reads 



dP 
'dt 



d'^P N -I dP 
+ 



r dr J ~^ \N / dx„ ~^ 



dP 



'l + 2/N)r'^-N P (48) 
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We define R = N ^/^r and get 



dP ed'P e{N-l) dP n : rr/v + 9^ ^^IP (AO) 

In the limit ^ 1 and e <^ 1, we drop the terms of 0(e), 0{N^^) and 
0{e/N), to obtain 



BP BP 

-g^ + Ril-R')g^ = [iN + 2)R'~N]P 



The method of characteristics now yields 

• 1 _ /?2 . jv/2 . 1 ^ l+Ar/2 
P(i?,t) ^ ' ^ 



G[e^\l-R^)/R' 



(50) 



(51) 



; Vi-i?2, 

Once again, this has to match smoothly to the small time solution in order 
to determine the unknown function G. Turning to Eq. ()48|) . the short-time 
solution is given by 



P(r < l,t < tc 



N/2 



L27re(e2* - 1)J 
1 

L27refe2* - 1)J 



g 2£{e2t-l) 



Ar/2 _ ivfl2 

g 2£{e2t-l) 



The matching is done by the choice 



(52) 



G[e^\l-R^)/R' 



R^ 



N 



1-P2 27re(e- 



2t 



N/2 



exp 



R^ 



N 



1 - i?2 2e(e2* - 1) 
(53) 

where exp(2t) has been replaced by (exp(2t) — 1) since exp(2t) ^ 1 and 
this replacement allows the matching with Eq. ()52p to be implemented. This 
leads to the answer 



N 

1 - i?2 2e(e2* - 1) 

(54) 

The discussion above is valid for a system of few degrees of freedom and 
as such the principle area of application is laser physics, where the laser oper- 
ates as a pump parameter crosses a threshold value and then the exponential 
growth is checked by a cubic nonlinearity in the Langevin equation. The vari- 
able is the electric field which, for circularly polarised light, can be taken to 
be two-dimensional. The experimental measurements of the time-dependent 
intensity clearly shows the existence of Suzuki regime. 



P{R,t) 



1 



J V/ z 



27re e^* - 1 



i?2 



exp 
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Of greater interest is the exploration of Suzuki scaling in the case of a 
field (Kawasaki et al., 1978; Bray, 1994), ie., a function of space and time, 
whose dynamics can be described by the Langevin equation 



d(j){f, t) 



6F 



+ f{r,t) 



(55) 



dt 5<i){r, t) 

where F is the "free energy" which can be written as (in a D-dimensional 
space) . 



(56) 



In principle, (j) may be a iV-component vector, and (j? = J2i=i 0? 

The growth of order that we have been considering crresponds to the 
parameter a being positive, if we work with the Fourier components (j){k, t) 
of 0(r, t), then 



d(j){k,t) 
dt 



a— 



e)(j){k,t)-X J2 <l){ki,t)(j){k2,t)<p{k-h-h,t)+f{k,t) (57) 



The probability distribution P{(j){k),t) satisfies 



dt 



E 



d 



P{-a + {k'y)(j){k',t) 



- d<j>{k'[ 

+\P J2 (P{h,t)(^{k2, t)(p{k' - ki - fcs, t) 
d^P 



d(f){k')d(f){-k' 



(58) 



The growth occurs for those fc- values which are smaller than a^^^, and 
our interest is in those wavenumbers alone. In the initial stages, where (f){x) 
is centred about = and is small in magnitude, we can drop the cubic 
term in the above equation and the time-development will occur according 

to P{(t>{k),t) ~ rjexp ^- I (t>{k') |2 /ee^''*^, where r = a - {k'f. In the 

intermediate time-zone, it is the e-term in Eq. ()58|l that needs to be dropped 
and we need to find the solution for P{(f){k),t) as we did in Eq. ljl^ . However, 
now it is more complicated. 

Simplification occurs if we look at the derivation of Eq.(0H) in slightly 
different manner. Returning to that case and setting e = 0, implies in the 
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Langevin picture solving the equation X = X — X^. The solution is X{t) = 
Ce'/{1 + C2e2*)V2 where C is a constant. If we call Xo{t) = Ce* = X(0)e*, 
then clearly X{t) = Xo{t)/[l + X^{t)]^/^ and X(0) = X{t)e-y[l - X^{t)Yl^. 
We thus have a prescription for going to X{Q) from X{t) and the evolution of 
P(X, 0) given as e~"^*^° will occur according to P{X,t) ~ \ 
In the present case, this requires the solution of 



dm t 

- (a — 



k^)(P{k)-\Y,<Pih,t)(P{k2,t)^ik-h-h,t) (59) 

ki,k2 

The solution (j){k) proceeds according to standard perturbative techniques. 
The simplification occurs in the large-time limit when we can identify a lead- 
ing term at every order. This allows the summation of the perturbation series 
to yield 

0(f,t) = 0o(r,t)/[l + 0^(f,t)]i/2 (60) 

with 

0o(r,t) = exp(t(a + V2))0(r-) (61) 
The inverse transformation is 

0(f) = exp(-t(a + V2))0(f)/[1 - 02(f)]'/' (62) 

In the absence of the diffusion term the probability follows a Liouville 
equation which is a conservation law for the probability. Hence, if the initial 
distribution is known, then the distribution at any time can be obtained by 
the inverse transformation shown in Eq.()62|). 

The evolution of this probability distribution can be pictured as follows. 
In the initial period, fluctuations everywhere start growing rapidly and at the 
same time diffuse over the distance (2rt)^/2, within which scale the fluctua- 
tions are strongly correlatad. Saturation sets in as the value of 02(f) at any 
point approaches unity. After a while, in the language of magnetism, the sys- 
tem breaks up into domains of size (2rt)^/2 with the saturation magnetisation 
at+1 or —1. 

Finally, we would like to point out connection with a field theory by 
writing the measure analogus to Eq.(jSII) as 

d/x ~ exp(-^[0]) Yl 1^<Paix, t) (63) 

where 

/oo 
ndt 
-co 

dt I d^x + v{4>. 



CO 



oo 



4e 



(64) 
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with the potential V{(j)a) coming from 

1 6U 6U 1 6^U 



V\ 



4e S(j)oi ^4>a 2 5(j)a5(j)a 



(65) 



With t = ir, the action is that for a quantum field theory with the action 



Sq, where 



drd^x 



4e dr dr 



+ v{4>o 



(66) 



lfU= (l/2)m^0a0a+(l/2)(V0Q,)^, ie., the Langevin potential is quadretic, 
Eq. ljU^ . prescribes 



Apart from constants, 



^ _ 2^2^ _d_d(pa ^ 
2 °'dxpdxji 



dxndx 



puxp 



(67) 



Sn 



drd^x 



1 



4e dr dr 2e dxp dxp 



4e 



dxftdx 



(68) 



which is a quadratic action and can be easily handled. 

It is when U is non trivial that one can generate non trivial Sq and what 
would be interesting is to consider a nontrivial S'g(e.g., the one corresponding 
to the decay of the false vacuum (Coleman, 1977)) and see if the corresponding 
Langevin dynamics can shed light on the quantum problem. 



4 A Class of Time-Dependent Potentials 

In this section, we will deal with potentials which are time-dependent but 
allow for the establishment of a final equilibrium state. This will be defferent 
from the time- dependences which are generally studied and have led to a 
wealth of interesting phenomena. These include the cases of diffusion over a 
barrier in the presence of harmonic force and diffusion over a fluctuating bar- 
rier. The hallmark of the former situation is the phenomenon of stochastic 
resonance (Benzi et al., 1981; Biittiker and Landauer, 1982; McNamara et al., 
1988), where the signal-to-noise ratio of the response to an applied force dis- 
plays a local maximum as a function of frequency. In the case of fluctuating 
barriers (Doering and Gadoua, 1992; Maddox, 1992; Ziircher and Doering, 
1993; Pechukas and Hanggi, 1994), the discovery that the mean first passage 
time has a minimum as a function of the correlation time characterising the 
fluctuation has prompted a wide variety of investigations. What we would 
like to present here is a toy model for yet another kind of phenomenon-the 
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global optimisation principle (Doye et al., 1999; Hunjan and Ramaswamy, 
2002) on an evolving potential energy landscape. In such problems one is in- 
terested in finding the minima of a multidimensional potential energy surface 
which constitutes the energy landscape in problems such as protein folding 
or the finding the ground state configuration of atomic or molecular clus- 
ters. An interesting observation in this context is that of Hunjan et al.(to be 
published), who have shown that a continuously and adiabatically varying 
potential assists approach to desired configuration at t — > oo by avoiding 
trapping in local minima. 

We want the final {t = oo) potential to be V2{x). We start from a different 
function and consider the time-dependent potential 

V{x, t) = V2{x) + b[Vi{x) - V2{x)] exp(-At) (69) 

which has the form b[Vi{x) — V2{x)] at t = and evolves to V2{x) at t = oo. 
For b = 1, Vi{x) is the initial shape which evolves to V2(x). Our goal is 
to study the approach to equilibrium in such a system. It is clear that 
as t ^ oo and V ^ V2, there is an equilibrium probability distribution 
Peg = exp(— V2/e) ai t = 00 for the Fokker-Planck equation. 

To study the onset of equilibrium when V is of the form shown in Eq. ()69|l . 
we make the usual substitution 

P{x,t) = 0(x,t)e-^("'*)/2^ 

(70) 

where 

Ay = Vi(x)-r2(x) (71) 

leading to 

^ = Ho4> + H,{t)^ (72) 



with 



and 



^o0 = e|^ + V,"0-l(l-,')^0 (73) 



Hi = b 



2 2e 2e 



4e 



In the above equation prime denotes differentiation with respect to x. 
The eigenvalues of Hq are non-positive and we can write 

Holpn = -Enlpn, with En > (75) 
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The ground state ipo has a space-independent part Vo{t). We will separate 
out this part and write the general solution of Eq. (f72|) as 

0(x,t) = 5^c.(t)e-^"*+/'^°(*')'^*V. (76) 
The usual techniques of Dirac's time dependent perturbation theory lead 

to 

Cn(t) = Cm{t) <m\H2\n> exp[~{Era - E„)t] (77) 

m 

where H2 = Hi — Vq. Perturbation theory proceeds by expanding 

Cn{t) = Cno{t) + bCnlit) + b^Cnl{t) + (78) 

Introducing it in Eq.(f77j) yields 

c;o = (79) 

cki{t) = Ycn,o{t) <m \ H \n> exp[-{E^ - En)t] (80) 

m 

Cn2{t) = Ycmi{t) < m \ H \ n > exp[-(K, - En)t\ (81) 

m 

and so on. We see immediately that Cno = constant, independent of time. 
Consequently, they are determined by the state of the system at t = 0. 

We focus on the where the final and initial potentials have qualitatively 
different structures. Our V2 will be a double well potential, while our Vi is 
the usual single well potential. Approach to equilibrium in the double well 
potential is strongly delayed, as we have seen in sec. 2, by the Kramers' 
time which is a long time scale coming from the possibility of noise induced 
hopping from one minimum to another. Thus V2{x) = ~\ ~\~ while 
Vi{x) = x^.With this choice 

Ho<f> = -er+(^^-^y (82) 
- (A/8e)x'(6 - x^) + (l/4e)x'(3 - x^)\l - e^*) 

(83) 

The eigenvalues spectrum of Hq is characterized, as we have seen before, 
by a set of close doublets as its low lying states. The separation within the 
doublet is exponentially small, while that between two doublets is of 0{1). 
The ground state is Eq = and the first exited state is the ground state of the 
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super symmetric partner of — ^ and is exponentially small. The second 
exited state has an eigenvalue close to 2 and hence we can safely approximate 
the dynamics of the low lying states as that of a 2-atate system. If (pQ and 
(pi are the two states and and 5 the eigenvalues, then we can write 

(Pix, t) = Co{t)(Po{x) + Ci{t)(Pi{x) exp(-5i) (84) 

Note that is the Kramers' time in the problem. The dynamics of Cq and 
ci is governed by 

Co =< (po\H'\(Po> e-^*co(t)+ <(po\H'\(Pi> ci{t)e-^' e'^' (85) 

ci = Co < 01 I I 00 > e''*e-^* + ci < 0i | i/' | 0i > e"^* (86) 

Now H' is even and hence < 0i | i7' | 0o >= 0, which decouples Co and 
Ci and we can easily integrate the above equations. If we drop terms like 
exp(— 2Ai), we get 

c'l = ci < 01 I ^' I 01 > exp(-Ai) (87) 

where 

H' = (3/2)(l - x") + (2xV4e)(3 - x'){l - x") - (A/8e)x'(6 - x") (88) 

The dominant contributions to both < 0i | if' | 0i > and < 0o | i/' | 0o > 
comes from near a; ~ 1 . The small difference between the two matrix elements 
comes from the term (3/2)(l— a;^) in H' , which is maximum near x — Q and at 
that point 0o 9^ but 0i ~ 0. After a set of straightforward manipulations, 
we see that 

P{x,t) = P,,{x)[l + f{x)eM-^efft)] (89) 

where 

5^ff = 6 + (90) 

where o; is a small number of the prder of S. 

The approach to equilibrium is now through a modified Kramers' time 
which is obtained from ^e//- Clearly 6eff > S and hence the new equilibration 
time is going to be smaller. Thus, in the toy model, we see a faster approach 
to equilibrium , which was the desired goal. 



18 



References 



[1] Benzi, R., Sutera, A. and Vulpiani, A. (1981). The mechanism of 
stochastic resonance. J. Phys. A, 14, L453. 

[2] Bernstein, M. and Brown, L.S. (1984). Supersymmetry and the Bistable 
Fokker-Planck equation. Phys. Rev. Lett., 52, 1933. 

[3] Bhattacharjee, J.K. and Banerjee, K. (1989). Kramers time in bistable 
potentials. J. Phys. A, 22, L1141. 

[4] Bray, A.J. (1994). Theory of phase ordering kinetics. Adv. Phys., 43, 
357. 

[5] Buttiker, M. and Landauer, R. (1982). Traversal time for tunneling. 
Phys. Rev. Lett., 49, 1739. 

[6] Caroh, B., Caroli, C. and Roulet, B. (1979). Diffusion in a bistable 
potential: a systematic WKB treatment J. Stat. Phys., 21, 415. 

[7] Coleman, S. (1977). Fate of the false vacuum: semiclassical theory. Phys. 
Rev. D, 15, 2929. 

[8] Doering, C.R. and Gadoua, J.C. (1992). Resonant activation over a fluc- 
tuating barrier. Phys. Rev. Lett., 69, 2318. 

[9] Doye, J.P.K., Miller, M.A. and Wales, D.J. (1999). The double-funnel 
energy landscape of the 38-atom Lennard- Jones cluster. J. Chem. Phys., 
110, 6896. See also Deaven, D.M. and Ho, K.M. (1995). Molecular ge- 
ometry optimization with a genetic algorithm. Phys. Rev. Lett., 75, 
288. 

[10] Hunjan J.S. and Ramaswamy, R. (2002). Global optimization by adia- 
batic switching. Int. J. Mol. Sci., 3, 30. 

[11] Hunjan, J.S., Matharoo, G.S., Sarkar, S. and Ramaswamy, R. (to be 
published) . 

[12] Kawasaki, K., Yalabik, M.C. and Gunton, J.D. (1978). Growth of fluc- 
tuations in quenched time-dependent Ginzburg-Landau model systems. 
Phys. Rev. A, 17, 455. 

[13] Kramers, H.A. (1940). Brownian motion in a field of force and the dif- 
fusion model of chemical reactions Physica, 7, 284. 

[14] Maddox, J. (1992). Nature, 359, 771. 

[15] McNamara, B., Wiesenfeld, K. and Roy, R. (1988). Observation of 
stochastic resonance in a ring laser Phys. Rev. Lett., 60, 2626. 



19 



[16] Pechukas, P. and Hanggi, P. (1994). Rates of activated processes with 
fluctuating barriers. Phys. Rev. Lett., 73, 2772. 

[17] Risken, H. (1984). Fokker-Planck Equation. Springer- Verlag, Berlin. 

[18] Schneider, T., Zannetti, M. and Badii, R. (1985). Stochastic simulation 
of quantum systems and critical dynamics. Phys. Rev. B, 31, 2941. 

[19] Suzuki, M. (1977). Scaling theory of transient phenomena near the in- 
stability point. J. Stat. Phys., 16, 11. 

[20] Zurcher, U. and Doering, C.R. (1993). Thermally activated escape over 
fluctuating barriers. Phys. Rev. E., 47, 3862. 



20 



